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ABSTRACT 

We use galaxy and dark halo data from the public database for the Millennium Sim- 
ulation to study the growth of galaxies in the iDe Lucia et al.l {2006) model for galaxy 
formation. Previous work has shown this model to reproduce many aspects of the 
systematic properties and the clustering of real galaxies, both in the nearby universe 
and at high redshift. It assumes the stellar masses of galaxies to increase through 
three processes, major mergers, the accretion of smaller satellite systems, and star 
formation. We show the relative importance of these three modes to be a strong func- 
tion of stellar mass and of redshift. Galaxy growth through major mergers depends 
strongly on stellar mass, but only weakly on redshift. Except for massive systems, mi- 
nor mergers contribute more to galaxy growth than major mergers at all redshifts and 
at all stellar masses. For galaxies significantly less massive than the Milky Way, star 
formation dominates the growth at all epochs. For galaxies significantly more massive 
than the Milky Way, growth through mergers is the dominant process at all epochs. 
At a stellar mass of 6 x 10 10 M©, star formation dominates at z > 1 and mergers at 
later times. At every stellar mass, the growth rates through star formation increase 
rapidly with increasing redshift. Specific star formation rates are a decreasing function 
of stellar mass not only at z — but also at all higher redshifts. For comparison, we 
carry out a similar analysis of the growth of dark matter halos. In contrast to the 
galaxies, growth rates depend strongly on redshift, but only weakly on mass. They 
agree qualitatively with analytic predictions for halo growth. 

Key words: galaxies: merger rate - galaxies: mass accretion rate - galaxies: star 
formation - cosmology: dark matter - cosmology: large-scale structure 



1 INTRODUCTION 

Galaxy mergers play an important role in galaxy forma- 
tion and evolution. They add new gas and stars. They 
drive gas motions which feed starbursts and central super- 
massive black holes, and, for comparab ly massi v e syst ems, 
they entirely restructure both galaxies. iToomrel (|l976l ) was 
the first to stress that the abundance of tidally distorted 
spirals in the nearby universe suggests that "star piles" 
produced by past interactions might accou n t for the ma- 
jority of observed elliptical galaxies. IWhitel (1 19781 ) carried 
out the first dynamically consistent 3-dimensional simula- 
tions showing that mergers do indeed produce remnants 
with a structure similar to that of ellipticals, a conclusion 
which has been reinforced by in creasingly realistic simu- 
lations of pu r ely stellar systems ([Farouki fc Shapiro! Il982l: 
iBarnes! Il988l ; iNaab fc Burkertl |2003l ). Inclusion of the gas 
component showed that a substantial fraction of the inter- 
stell ar medium should be driven t o the centre in major merg- 
ers l|Negroponte fc White! Il983l ; H arncs fc Hernauisti ri991: 
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IMihos fc Hernauistlll99r3 ). This work supported the identi- 
ficati on of ultraluminous infrared galaxies as merging sys- 
tems l|Sanders et al.|[l988l ) but led to remnant galaxies with 
cores which are denser than observed ellipticals. Recent work 
suggests that this contradiction may be resolved by strong 
AGN- or starburst-generated winds which expel a large 
fraction of the gas fro m the galaxy (jSpringel et al.ll2005al ; 
|Pi Matteo et alJboOot ). Work on mergers of unequal galax- 
ies suggests that while such merg ers may not greatly alter 
the structure of th e larger system (jVelazquez fc Whitell 19991 : 
lAbadi et al.ll200&D they can nevertheless stimulate substan- 
tial rearrangements of its gas with associated star formation 
and AGN activity (iMihos fc Hernquist|[l99i ). 

In the standard ACDM cosmology structure forms hi- 
erarchically. Small dark matter halos form first and then ag- 
gregate into progressively larger systems. At any given time 
cosmic matter is distributed over nonlinear objects spanning 
many decades in mass, and growth is driven by merging 
with similar halos, by accretion of much smaller halos and 
of diffuse material, and by des truction by infall onto larger 
halos (e.g. lLacev fc Colelll993l ). The situation is made more 
complex by the fact that the inner cores of halos often sur- 
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vive as long-lived substru cture within the larger objects by 
which they are accre ted l|Ghigna et all 1 199a Im oore et all 
Il999l : iGao et al.ll2004l ). Gala xies form at th e centr es of ha- 
los in the way suggested by IWhite fc Reesl { 19781 ) and are 
swept along with the growth of dark matter structure. They 
gain stars through formation from their interstellar medium, 
which may be replenished by infall from their surroundings, 
and by incorporating the stars of galaxies which merge with 
them. The interaction between these processes drives the 
overall evolution of the population and cannot be followed 
without treating the associated baryonic astrophysics (gas 
condensation, formation and evolution of stars and black 
holes, feedback from supernovae and AGN, chemical enrich- 
ment, production of observable radiation etc.). 

Early studies of the evolution of the galaxy population 
embedded simplifed treatments of this baryonic physics in 
Monte Carlo realisations of the merger tre es associated with 
the formation of individual dark halos jK auffman n et all 
Il993l : ICole et all 1 1994 120001 : ISomerville fc Primackl 1 19991 ). 
The spatial distribution of galaxies could then be studied 
using the halo distri bution from an N-body simulation of 
struc ture formation (|Kauffmann et all Il997l : iBenson et all 
2000). Improvements on this scheme have used higher res- 
olution N-body simulations so that the merging trees can 
taken directly from the simulation itself, thereby allow- 
ing the evolution of the galaxy po pulation to be fo l lowed 
in a single consistent simulation (|Kau ffmann e t al. 
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Kang ct al. 2005). A parallel approach 



has followed the dynamics of diffuse gas (in particular, as- 
pects of the gas condensation and galactic wind processes) 
by adding a hydrodynamic scheme to the N-body treat- 
ment of dark matter while continuing to treat star formation 
and evolution by semi-analytic means (|Cen fc Ostriker 1992; 
iNavarro fc Whitell 19941 ; iKatz et al.lll99rj ). The development 
path here has involved continual improvement of the simu- 
lation schemes to increase resolution and to treat the acces- 
sible physics more realistically (e.g. | Cen fc Ostrikerl |2000| ; 
ISpringel fc Hern quist 200l; EE ommer et all 120061 ) . Recent 
work in both approaches has focussed on how feedback from 
AGN may re gulate the formation and evolution of their 
host galaxies (ISpringel et all l2005al; [Pi Matteo et all [20051 ; 



ICroton et aHl2006l ; iBower et alll200e 

This body of work has demonstrated that while galaxy 
mergers are an important aspect of the evolution of the 
galaxy population, t hey do not simply pa rallel the merg- 
ers of dark halos. As I White fc Reesl (| 19781 ) stressed, galax- 
ies must remain distinct after the merger of their hal os if 
we are to understand the formation of galaxy clusters. iFalll 
(1979) noted that late-type giant galaxies cannot have un- 
dergone recent major mergers since these would destroy 
their ste llar disk s. Whi le many more recent studies have 
followed iToomrd |l976) in arguing that massive ellipti- 
cal galax ies assembled relatively rec e ntly through merg- 
ers (e.g. iKauffmann fc Charlotl Il998l ; Ivan Dokkuml 120051 ; 
|Pe Lucia et al.ll2006l ) other authors have used the age and 
uniformity of their stellar populations and their apparently 
undiminished abundanc e at high redshift t o argue against 
such late assembly (e.g. ICimatti et all l2006). Observational 
estimates of merging rates, based primarily on counts of 
very close pairs of galaxies, or of morphological evidence 
for recent merging, have varied widely due to uncertainties 



in th e associated timescales (|Le Fevre et allfeood : ILin et all 
120041) . In addition, attempts to measure the evolution of the 
merger rate, usually parametrised as proportional to (l + z) Q 
ha ve obtained values for the exponent a ranging from to 
6. dBell et al.ll200rj; ICaxlberg et all^OOol; iPatton et al.ll2002l ; 
IConselice et al.ll2003l ; iBundv et al.ll2004l ; ILin et al.ll2004l ). 

In the present paper we analyse the build-up of 
the galaxy populat i on in the galaxy formation model of 
iDe Lucia fc BlaizotJ (|2007h which i s implemented on th e 
very large Millennium Simul ation (ISpringel et all l2005bl ) . 
This model updates that of (|Croton et al 20061 ) and is a 
reasonable match to the clustering and to many of the 
systematic properties of the local galaxy population. It is 
also consistent with most available data at high-redshift 
( Kitz bichler fc Whitell2007l ). For our purposes, this provides 
a physically consistent and observationally plausible imple- 
mentation of galaxy formation within the dynamical frame- 
work of ACDM. It can therefore be used to explore the dif- 
ferences between galaxy growth and dark halo growth in 
this structure formation model. We use the public database 
containing the properties of the dark halos and the galax- 
ie^] to construct mean growth rates for galaxies through 
major mergers, through minor mergers and through star- 
formation, each as a function of galaxy mass and of redshift, 
and we compare these with analogously defined growth rates 
for dark halos. 

Our paper is organized as follows. In Sec.[2]we introduce 
the Millennium Run and the prescriptions used to simulate 
galaxy formation using merger trees built from it. Sec. [3] 
presents our analysis of the mass and redshift dependence 
both of the major merger rate and of growth rates through 
major and minor mergers as well as through star formation. 
In Sec. |4j we discuss the corresponding properties of dark 
halos (defined here as FOF groups) and contrast them with 
our results for galaxies. Conclusions and discussions are pre- 
sented in Sec. [5] 



2 THE SIMULATION AND THE GALAXY 
FORMATION MODEL 

The galaxy catalogue used in this paper was produced us- 
ing a "hybrid " technique: a large N-body simulation was 
first carried out to define the evolution of the dark matter 
distribution, and then a suite of semi-analytic prescriptions 
was implemented in order to simulate the formation and 
evolution of galaxies within a stored "forest" of (sub) halo 
merging trees constructed from the original simulation. A 
detailed description of the Millennium Sim ulation and of 
the ga l axy formation mod e l can be found in (jSpringel et all 
l2005bl ; ICroton et ai1l2006l ; IDe Lucia fc Blaizotl l2007l ). Here 
we summarize the main simulation characteristics and the 
way the halo merger trees were constructed, as well as those 
aspects of the galaxy formation modelling that are relevant 
to our study of galaxy growth. 
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vals thereafter. At each time, the simulation code produced 
a friends-of-friends group catalogue on the fly by linking to- 
gether particles separa ted by less than .2 of the average in- 
terparticle separation (iDavis et alj 1985th S ubsequently, the 
SUBFIND algorithm (|Springel et al.ll200lh was used to di- 
vide each FOF group into a disjoint set of self-bound subha- 
los. These subhalos are the basis for the merger trees, which 
are defined by linking each subhalo from a given output time 
to one and only one descendent at the following output time. 
When studying the growth of dark halos in Sec. |4j we define 
a halo as an FOF group and we estimate its mass as the sum 
of the masses of all its subhalos. This typically loses a small 
amount of "diffuse" material which was bound to none of 
the subhalos. This is not significant for our purpose here, 
and this definition was convenient, since the original FOF 
halo mass was not stored in the (sub)halo database when we 
carried out this project. More importantly, this mass defi- 
nition allows us to deal in a straightforward way with the 
problem that simulated halos, unlike our simulated galaxies 
or the halos considered in simplified models for halo growth, 
not only merge but can also fragment. Mass from a single 
FOF halo can thus contribute to several FOF halos at some 
later time. 



1 +z 



Figure 1. The specific rate of formation of galaxies through 
major mergers as a function of redshift. The 7 curves refer to 
product galaxies with stellar mass in 7 disjoint ranges, identified 
by labels with units of 10 10 Mq. A galaxy is considered to have 
just formed through a major merger if its two largest progenitors 
in the preceding Millennium Run output differ by less than a 
factor of 3 in stellar mass. The dimensionless measure of merger 
rate used here is the fraction of all galaxies in the given stellar 
mass bin at redshift z which form through a major merger per 
unit time, multiplied by the age of the Universe at redshift z. 
Error bars give Poisson uncertainties derived from the number of 
new merger products in each bin at each redshift. The probability 
that a galaxy has just formed through a major merger is a strong 
function of stellar mass, but a weak function of redshift. 



2.1 The simulation 



The Millennium Simulation (|Springel et alj|2005bh used in 
this study is the largest simulation of cosmic structure for- 
mation so far carried out. It adopts the concordance ACDM 
cosmology and follows N = 2160 3 particles from redshift 
z = 127 to z — in a comoving box of side-length bOOh" 1 
Mpc. This volume is large enough to investigate rare ob- 
jects such as quasars and rich clusters of galaxies, yet, with 
N — 2160 3 particles, has a dark matter particle mass of 
only 8.6 * 1O 8 M0, allowing the galaxy formation model to 
follow the formation of all galaxies more massive than the 
Small Magellanic Cloud. The assumed cosmological param- 
eters are Q m = 0.25, fi b = 0.045, h = 0.73, ft A = 0.75, 
n = 1, and erg = 0.9, where the Hubble constant is parame- 
terized as usual as Ho — WOhkms' 1 Mpc" 1 . These param- 
eters are consistent with a combined analysis of the 2dF- 
GRS llColless et alJ I2OO1T ) and the first-year WMAP data 
l|Spergel et aUl2003 ). 

During the simulation, the full particle data were stored 
at 64 output times approximately logarithmically spaced 
from z = 20 until z = 1 and at approximately 200 Myr inter- 



2.2 Merger rates 

In the galaxy formation models implemented on the Millen- 
nium Simulation, a galaxy begins to condense at the centre 
of a halo as soon as it is identified as a persistent object with 
more than 20 dark matter particles. As the halo grows, so 
does the galaxy at its centre, forming stars at a rate gov- 
erned by its cold gas content and by empirically determined 
star formation "laws" . The halo may merge into a larger sys- 
tem, becoming an independent subhalo orbiting within the 
FOF group. The galaxy is now considered a satellite, losing 
its supply of fresh gas, and perhaps ceasing to form stars if 
it uses up its available interstellar medium. Dynamical fric- 
tion effects bring the orbit of the subhalo ever closer to the 
centre of its parent, and tidal effects strip away its outer 
regions until eventually it may be disrupted entirely (or at 
least drop below the resolution limit of the simulation). At 
this point the galaxy is associated to the most-bound par- 
ticle of the subhalo at the last time it was identified and is 
marked as a candidate for merging with the central galaxy 
of the parent halo. The merger takes place one estimated 
dynamical friction time later. 

Galaxy mergers may trigger strong star formatio n. In 
the galaxy formati o n mo del of ICroton et al.l (|2006t ) and 
|Pe Lucia fc BlaizotJ d200Tt> which we ana lyse here, a recipe 
similar to that of ISomerville et ahl (|200lT ) is adopted to de- 
scribe starbursts during minor mergers. In this model, a frac- 
tion eburst of the cold gas of final galaxy is converted into 
stars, where 



eburst = 0.56 * ( 



M. 



satellite \Q.7 



M r , 



A major merger is assumed to occur whenever the two galax- 
ies differ by a factor of less than 3 in baryonic mass. In such 
a merger the starburst is assumed to convert a large frac- 
tion of the cold gas into stars and to eject the rest from the 
galaxy. The remnant of such a merger is assumed to be an 
elliptical galaxy. It may, however, grow a new disk if gas is 
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able to cool from the surrounding halo, and in this case the 
merged system becomes the bulge of a larger spiral galaxy. 

In this galaxy formation model, central galaxies are 
treated differently than satellites. Only central galaxies are 
fed new material by cooling from the hot atmosphere of their 
halo, by direct infall of cold gas, or by merging of satel- 
lites. No new material accretes onto satellite galaxies, so 
that their star formation terminates when their cold gas is 
used up. Gas accretion processes depend strongly on time 
and on galaxy mass. At early times and in low-mass galax- 
ies gas cools substantially more efficiently than in high-mass 
systems and at late ti mes. In additio n , an i mportant inno- 
vation in the model of ICroton et al. (2006) (and included 
here) is a treatment of "radio mode" feedback. This assumes 
that if the central galaxy has a supermassive black hole and 
sits at the centre of a static hot atmosphere, then radio ac- 
tivity will prevent further cooling of hot gas. This resolves 
the long-standing "cooling flow problem" and ensures that a 
massive elliptical at the centre of a group or cluster does not 
grow a new disk and so remains "red and dead" . As a result 
the only significant growth mode for high mass galaxies is 
through merging. 

In this study, we consider all galaxies in the Mil- 
lennium/DeLucia database with stellar mass between 5 x 
10 9 M Q and 6.4 x 1O 11 M . Although the galaxy catalogues 
are nearly complete to a mass at least 5 times lower than 
this, we want to be able to resolve the recent merging his- 
tory of each system and so we adopt this more conservative 
limit. This choice leaves us with a total 81896686 galaxies 
(summed over all redshifts). To investigate the mass depen- 
dence of galaxy growth, we divide this sample into seven 
mass bins, each a factor of 2 wide. The highest mass bin 
contains the smallest number of galaxies, a total of 22827 
systems. 



3 GALAXY GROWTH RATES 

Growth in the stellar mass of galaxies occurs through two 
processes: conversion of gas into stars (either quiescently or 
in a starburst) and the addition of stars through mergers. In 
this section we mine the publicly available database to study 
the interplay between these processes. We begin by study- 
ing how the rate of major mergers depends on the mass of 
the product galaxy and on redshift. We then compare mean 
galaxy growth rates due to this process to mean growth rates 
due to all mergers (major and minor) and to star formation. 
For each galaxy in the database at each time, we define the 
main progenitor at the previous stored time to be the pro- 
genitor with the largest stellar mass. If a galaxy has more 
than one progenitor at the earlier time, then it has under- 
gone a merger between the two times. If m of the other 
progenitors differ from the main progenitor by less than a 
factor of 3 in stellar mass, then the galaxy is assumed to 
have had m major mergers in this time interval. 

We define a dimensionless major merger rate per galaxy 
as a function of redshift and stellar mass through 



10.00 



R(M*,z) = 



N major {M t ,z)/5t{z) 
N gat (M,,z)/t(z) 



(1) 
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where N ga i (M* , z) is the number of galaxies in the simula- 
tion at redshift z and with stellar mass in a chosen inter- 



Figure 2. The relation between the stellar mass of galaxies and 
their specific formation rate through major mergers. The rates 
given here average the data plotted in Figure 1 over the redshift 
range z sC 2 (except for the highest stellar mass bins where 
there are insufficient objects to determine a rate at the higher 
redshifts). The error bars indicate the rras fluctuation in rate 
over these redshift intervals. Clearly, the probability that a galaxy 
has just formed through a major merger increases approximately 
linearly with stellar mass in this galaxy formation simulation. 



val centred on M*, N merger (M*, z) is the number of these 
galaxies which have had a major merger since the last stored 
redshift z p (z) (a galaxy which has had m major mergers is 
counted m times), 5t(z) is the time interval between z p and 
z and t(z) is the age of the universe at z. Hence R(M*,z) 
is the fraction of galaxies of stellar mass M* formed per 
Hubble time through major mergers at redshift z. 

Fig.[TJshows major merger rates estimated in this way as 
a function of redshift for seven intervals of stellar mass, each 
a factor of 2 wide. We plot Poisson errors on our estimates 
which are determined entirely by the number of merger rem- 
nants N m ajor(M* , z) found at each time. At low redshift 
(z < 2) our dimensionless rate depends remarkably weakly 
on redshift. For large stellar mass, any variation is within 
the noise. For smaller M* there is a gentle rise (by less than 
a factor of 2) from low redshift out to z ~ 1.5, followed by a 
decline at higher redshift. On the other hand the dependence 
of R(M t , z) on stellar mass is very strong. The probability of 
formation through major mergers is about 40 times higher 
for the most massive galaxies we consider than for the least 
massive galaxies. Galaxies comparable in mass to our Milky 
Way form through major mergers at a rate of about 25% of 
the population per Hubble time, while for galaxies with a 
stellar mass ~4x 10 11 M Q the corresponding rate is about 
8 times higher. 

To see more clearly the stellar mass dependence of the 
specific rate of formation through major mergers, we plot in 
Fig.[2]the relation between {R(M*, z)) and stellar mass. Here 
we have averaged the data of Fig. [TJover the redshift interval 
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Figure 3. Dimensionless mean growth rates for galaxies as a function of redshift for the 7 different stellar mass bins of Figure 1 and 
for different growth modes. For each galaxy at each output time, the fraction of its stellar mass gained in a particular mode since the 
previous output is divided by the time between outputs and multiplied by the current age of the Universe. The result is then averaged 
over all galaxies in the chosen mass bin and plotted against output redshift. The different curves represent stellar mass growth through 
major mergers (orange) through all mergers (green) and through star formation (black). The stellar mass ranges in the labels for each 
panel are given in units of 10 10 Mq. 
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from z = 2 to (or, for massive galaxies, over redshifts where 
there are more than 15 major mergers in total). Error bars 
show the rms variation in the rate over the redshift range 
used. The relative formation rate through major mergers 
is approximately proportional to stellar mass R oc M* , al- 
though the plot suggests a more complex behaviour with an 
initial steepening towards higher mass followed by a (possi- 
ble) saturation at the highest mass. 

Galaxies grow not only through major mergers, but also 
through minor mergers and through star formation. In or- 
der to compare the relative importance of these processes, 
we now calculate mean growth rates for galaxies in each of 
these channels as a function of stellar mass and redshift. 
In analogy to equation (1) we define mean dimensionless 
growth rates due to major mergers, to all mergers and to 
star formation as 



Rn 



M m ajor{M„,z)/St(z) 

M gal (M*,z)/t(z) 



R 



m,merger 



,(M t ,z)/St(z) 



R m ,gas(M t ,z) = 



M ga i(M*,z)/t(z) 
M gas (M,,z)/St{z) 



M ga i(M*,z)/t(z) 



(2) 



(3) 



(4) 



where M ga i (M* , z) is the total stellar mass of all galaxies at 
redshift z with individual stellar masses in the bin centred on 
M„, and M with subscripts 'major', 'merger' and 'gas' indi- 
cates the total stellar mass added to the main progenitors of 
these galaxies since the previous output time through major 
mergers, all mergers and star formation, respectively. This 
includes star formation over this time interval in all the pro- 
genitor galaxies, as well as in quiescent and in merger-related 
starburst modes. St(z) and t(z) have the same meaning as 
before. These rates represent the recent growth of galaxies 
prior to the time they are observed in terms of the fractional 
increase in their stellar mass per current Hubble time occur- 
ring in each of the three modes. For example, R m , g as > 1 
represents a class of galaxies whose recent average star for- 
mation rate exceeds their past average star formation rate. 

In Fig. we plot these growth rates as a function of 
redshift for the same 7 bins of stellar mass already illustrated 
in Fig. [1] The orange curves give the dimensionless growth 
rate through major mergers and so are very similar to the 
curves already plotted in Fig. [1] Indeed, the ratio of the two 
is just the average of the ratio of the stellar mass of the 
smaller galaxy in a major merger to the stellar mass of the 
merger product. Thus, the dimensionless growth of galaxies 
through major mergers also depends little on redshift but 
strongly on stellar mass (as in Fig. Only for the most 
massive galaxies does R ma jor approach unity; for galaxies 
of Milky Way mass it is around 10% at all redshifts. 

The green curves in Fig. [3] give mean growth rates due 
to all mergers. For all but the more massive galaxies at the 
lowest redshifts, these curves lie more than a factor of 2 
above the major merger curves. The difference between the 
two curves increases with increasing redshift in all cases. 
Thus, minor mergers are generally more important for in- 
creasing the stellar mass of galaxies than are major merg- 
ers. For small mass galaxies at high redshift the ratio of the 
two growth rates can be an order of magnitude. For galaxies 
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Figure 4. For comparison with Fig. [T] we plot redshift against 
the specific formation rate of FOF groups through major mergers, 
averaged over groups in 6 different mass bins as indicated by label 
color. The mass unit here is 10 10 -Mq. 



with masses above 10 n Mo (in the model this represents the 
classical giant elliptical population) merging dominates the 
growth rates at redshifts z < 2, and major mergers account 
for more than half of the total stellar mass growth at low 
redshifts. In the highest stellar mass bin the relative impor- 
tance of major and minor mergers is slightly different; these 
objects are the Brightest Cluster Galaxies investigated in 
detail by De Lucia & Blaizot (2007). 

Finally, the black curves in Fig. [3] give mean growth 
rates due to star formation as a function of redshift. These 
are constructed by averaging all the star formation between 
two output times in all the progenitors of the galaxies in 
each mass bin. As a result, they include quiescent star forma- 
tion both in the main galaxies and in smaller galaxies which 
merge with them, as well merger-induced starbursts. Unlike 
the growth rates due to mergers, they increase monotoni- 
cally and relatively steeply towards high redshift, roughly as 
one power of (1 + z) on average, although the slope decreases 
with redshift at low stellar mass and increases with redshift 
at high stellar mass. At the present day Rm,gas(M*, 0) is a 
decreasing function of M* and is always below unity. Thus, 
galaxies of all stellar masses are, on average, currently form- 
ing stars at less than their past average rate. For galaxies 
of Milky Way mass, the mean star formation rate at z = 
is about 15% of the past average; this ratio drops to very 
small values for more massive systems. 

This behaviour is well known in the real Universe and is 
often taken as evidence for "downsizing"; massive galaxies 
seem to complete most of their star formation at higher red- 
shift than low mass systems. Somewhat surprisingly, how- 
ever, this ranking of dimensionless growth rate holds at all 
redshifts, not just at z = 0. In this model there is no red- 
shift at which high stellar mass galaxies are growing faster 
(in relative terms) than less massive systems. Except for the 
highest mass bin (where galaxies form almost exclusively 
through multiple mergers) the dimensionless growth rates 
due to star formation exceed unity at sufficiently high red- 
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shift for galaxies of all stellar mass. This remains true to 
lower redshift for lower stellar mass. 

If we compare the mean growth rates due to star for- 
mation with those due to (all) mergers, we see that, except 
at the highest stellar masses, star formation dominates at 
sufficiently high redshift. This is true all the way down to 
z — for galaxies less massive than the Milky Way, but 
for higher mass systems mergers are the dominant growth 
mode at low redshift. It is interesting that the Milky Way 
mass, which is also approximately the characteristic stellar 
mass at the knee of the galaxy luminosity function, marks 
the boundary between the two regimes. This is not a coinci- 
dence. It is built into the model by the physical assumptions 
required to get a good fit to the observed galaxy luminos- 
ity function. In low-mass systems cooling is very efficient 
and supernova feedback has to be invoked to prevent over- 
production of stars. Even with such feedback, a significant 
fraction of the baryonic material gained by small halos is 
turned into stars, and most of this accreted material is asso- 
ciated with objects which were too small to contain stars of 
their own. Hence star formation is a more effective growth 
mode than merging. At Milky Way mass, cooling is still ef- 
ficient, particularly at early times, and supernova feedback 
is less effective in preventing star formation. On the other 
hand, much of the infalling material is in objects which are 
massive enough to contain substantial numbers of their own 
stars. Thus stellar mergers become competitive with star 
formation. For higher stellar masses, the model invokes "ra- 
dio mode" AGN feedback to suppress cooling and star for- 
mation. The steep quasi-exponential tail of the stellar mass 
function is then populated almost exclusively by mergers. 



4 GROWTH RATES FOR FOF GROUPS 

As discussed in Sec. [TJ the relation between galaxy mergers 
and mergers of their host halos is less straightforward than 
one might expect. In this section we investigate merger and 
growth rates for dark halos in a way which allows direct 
comparison with the results presented for galaxies above. 
For the purposes of our study it is convenient to identify 
dark halos as the friends-of-friends (FOF) groups initially 
identified in the Millennium Simulation, and to approximate 
the mass of each FOF group by the sum of the masses of its 
identified subhalos. This loses the mass of a certain number 
of "diffuse" particles which are not bound to any subhalo, 
but this systematic is relatively small for most halos and is 
of no consequence for our analysis. This scheme provides a 
straightforward way for us to deal with the problem that 
simulated halos, unlike those in extended Press-Schechter 
models (e.g. Lacey & Cole 1993) or the galaxies discussed 
above, often fragment into pieces which can become parts of 
different halos at a later time. This means that the progeni- 
tors of an FOF halo may include only part of an earlier FOF 
halo. Tracking individual subhalos allows us to account for 
this, since the Millennium halo database is set up so that 
each subhalo has a unique descendent, ensuring that the 
progenitors of an FOF group are a unique set of subhalos 
which may form all or part of several FOF halos. 

We bin the FOF groups at each redshift according to 
mass, with each bin spanning a factor of 3.8 in mass. The 
six bins for which we present results then correspond very 



roughly to halos whose central galaxies lie in the upper six 
stellar mass bins of figures [TJ and [3] 

In Fig 3J we plot the redshift dependence of the specific 
rate of formation of FOF halos through major mergers for 
our 6 bins of halo mass. A halo is defined to have just un- 
dergone m major mergers if its progenitor subhalos at the 
previous output come from at least m + 1 different FOF 
halos, and if the total subhalo mass coming from m of the 
subdominant FOF progenitors is more than a third of that 
coming from the main FOF progenitor. This merger count 
can then be used to define a merger rate in direct analogy 
to equation (1). Figure 3] can be compared directly with the 
6 most massive cases of Fig. [TJ The behaviour is quite dif- 
ferent, however. In Fig. [4] there is a strong and monotonic 
dependence of formation rate on redshift, but there is lit- 
tle dependence on halo mass. This is the exact contrary of 
what we found for galaxies, where the mass dependence was 
strong and the redshift dependence weak. The redshift de- 
pendence of these curves is reasonably well described as a 
simple proportionality to (1 + z). For all masses the rates 
exceed unity for all but the lowest redshifts. Recall that in 
Fig. [TJ we found the coresponding rates for galaxies to exceed 
unity only for the most massive systems. Major mergers are 
thus a much more significant growth mode for dark halos 
than they are for most galaxies. 

Fig [5] shows dimensionless growth rates for FOF halos 
as a function of redshift for the same 6 halo mass bins. These 
rates are defined in exact analogy to equations (4) through 
(2) and refer to growth through major mergers (orange), 
through all resolved mergers (green), and through accretion 
of "diffuse" particles (i.e. simulation particles not assigned to 
any FOF halo with more than 20 particles; the black curve). 
Again the growth rate through major mergers parallels the 
specific formation rate already plotted in Fig. [3] the ratio 
of the two is just the average mass of the smaller partner in 
a major merger in units of the final halo mass. Both the 
growth rate through major mergers and the growth rate 
through all (resolved) mergers are near power-laws of similar 
slope. The growth rate through all resolved mergers exceeds 
that through major mergers by a larger factor for high-mass 
halos than for low-mass ones. This primarily reflects the fact 
that the resolution limit of the simulation corresponds to a 
much lower mass ratio limit for identifying a merger in the 
former case. This is not the whole story, however, as one 
can see by the fact that the diffuse accretion rate depends 
differently on redshift than the other growth rates. Hence the 
growth of objects of given mass is more strongly affected by 
accretion of diffuse material at early times than at late times. 
In addition, comparing the major merger growth rates (the 
orange curves) with the total growth rates (the sum of the 
green and black curves), one sees that while at high masses 
and at early redshifts major mergers account for about 15% 
of the total growth rate, for small objects at late times they 
account for a larger fraction of the growth. Note that at all 
redshifts and for all masses, accretion of "diffuse" particles 
accounts for at least 30% of the total growth. 



5 SUMMARY AND DISCUSSION 

We have used public ly available data for the 
|Pe Lucia fc Blaizoj (|2007l ) model of galaxy formation 
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Figure 5. For comparison with Fig. [3] we plot dimensionless mass accretion rates for FOF groups as a function of redshift for the same 
6 mass bins as in Fig. [4] The different curves give the mean mass accretion rate due to major mergers (orange), to all mergers (green) 
and to accretion of diffuse particles (black). The mass unit for the labels in each panel is 10 10 Mq. 



to study the relative importance of merging and of star 
formation for the growth of galaxies. This model is based on 
stored halo merging trees for the Millennium Simulation, 
a very large simulation of the evolution of the dark matter 
distribution in a ACDM cosmology. It is consistent with a 
wide variety of observational data on the properties and 
clustering of galaxies both at low and at high redshift. 
We thus expect its behaviour to give at least a qualitative 
indication of the balance needed between the various modes 
of galaxy assembly in any successful model in the ACDM 
context. A particular goal of our study has been to contrast 
the roles of merging in galaxy and dark halo evolution. 

The most striking result from our study is that forma- 
tion through merging depends in a completely different way 
on mass and redshift for our two classes of object. Recent 
formation through a major merger is almost equally likely 
for halos of all masses at any given time, but is substantially 
more likely at early times than it is today. For galaxies, on 
the other hand, the likelihood of recent formation through 



a major merger is a strong function of stellar mass, but de- 
pends at most weakly on redshift. In addition, halos of all 
masses have grown more rapidly through mergers than all 
but the most massive galaxies. A little reflection shows that 
these differences are required by the facts that a galaxy clus- 
ter is considered as a single dark matter halo but contains 
many distinct galaxies, and that the stellar mass function 
for cluster galaxies differs little from that of the Universe 
as a whole. This implies that the build-up of massive halos 
through mergers cannot be paralleled by merging of the as- 
sociated galaxies. Merging plays a much less important role 
(though still significant) in galaxy growth than in dark halo 
growth. The high rates of recent merging found for the most 
massive galaxies are a selection effect. Only through merg- 
ing can galaxies attain such high masses. This is also the 
reason why the most massive galaxies are usually ellipticals, 

A second striking result from our study is the increas- 
ing importance of star formation with increasing redshift for 
galaxies of all masses. At low redshift we find the observed 
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result that mean specific growth rates through star forma- 
tion are smaller in high-mass galaxies than in low-mass ones, 
but it turns out that this result also h olds a t high redshift. 
According to the I5e~ Lucia fc Blaizotl (120071 ) model there is 
no redshift where the specific star formation rate of massive 
galaxies significantly exceeds that of low mass systems. Indi- 
vidual objects may be experiencing dramatic starbursts, but 
averaged over the population of all objects of given stellar 
mass, the prediction is that the mean specific growth rate 
through star formation is always a decreasing function of 
stellar mass. 

Only at redshifts below one and for galaxies compara- 
ble to or more massive than the Milky Way does the growth 
rate through mergers exceed that through star formation. 
This corresponds nicely to the "transition stellar mass" at 
which the stellar populations and the structural parame- 
ters of local galaxies switch from being predominantly star- 
forming and disk- l ike to predominantly old and spheroidal 
l|Kauffmann et al.l l2003 l. This agreement is, of course, in 
part a consequence of the tuning of the parameters of the 
galaxy formation model to fit observation. 

A less surprising but still interesting result is that 
merger-related growth for objects of all stellar masses and at 
most times is roughly equally divided between what we have 
defined as major and minor mergers. Clearly our separation 
at a progenitor stellar mass ratio of 3 to 1 is arbitrary. If we 
had chosen 5 to 1, major mergers would have dominated in 
most cases. If we had chosen 50 to 1, minor mergers would 
have been unimportant. Clearly the accretion of the LMC 
will make a much more significant change to the Milky Way's 
stellar mass than the addition of all the Dwarf Spheroidals, 
and this in turn will be dwarfed by the impending merger 
with M31! 

As we now show, the FOF halo behaviour we find is 
at least qualitatively consist ent with the predictions of EPS 
theory l|Lacev fc ColeHl993T ). The analytical expression for 
the probability that a mass element which is part of a halo 
of given mass M2 at time fa is part of a halos of (smaller) 
mass Mi at the earlier time ti is 

f(Sl,LU 1 \S2,LU2)dSl 

_ UJ! - L0 2 r (qji - ^2) 2 1-JO 

~ (2^)V2( 5l _ S 2f/ 2 6XP[ 2{Sx - S 2 y 1 

where Si. 2 is the rms linear density fluctuation (extrapo- 
lated to z = 0) in spheres containing a mean mass Mi, 2, 
wi,2 = 8 c o/ D(zi : 2) is the redshift- dependent critical density 
for collapse, D(z) is the growth factor of linear fluctuations, 
and S c o ~ 1.69 is a constant. By taking the limit as £2 tends 
to t\ (so u)2 — wi tends to 0) and integrating over Si , we can 
get the dimensionless merger rate per product halo: 

P(M high , M [ow ,UJi\S2,0J2) 

dt/t 

1 r S{Mh ^ h) m 2 t*dw/dt 

~ 2 's { M lom) Mi( 27 r)V 2 (Si-S 2 )3/ 2 1 

Setting M tow = ^ and M high = ^j 2 - we get the major 
merger rate which is seen to evolve with time as t * duo/dt. 
In a Einstein de Sitter universe D(z) oc (1 + z) -1 oc t 2/3 and 
thus t * dui/dt oc (1 + z), roughly reproducing the behaviour 
we get for the major merger rate of FOF halos in the Mil- 



lennium Simulation. In the ACDM cosmology, the formula 
is more comp l ex bu t is quantitatively similar. As shown by 
ICarroll et ail (|l992h . D(z) = g(z)/[g(0)(l + z)] where 

g(z) » 5/2n m [Qi( 7 _ n A + (1 + On/2)(l + ^a/70)]" 1 

and Q, m (Qa) is the density parameter of matter (dark en- 
ergy). We plot t * dui/dt against 1 + z for the two cases in 
Fig 6 to illustrate the size of the expected differences. 

The same formalism also allows the dimensionless mass 
accretion rates through mergers and/or smooth accretion to 
be expressed as 

Pm (M h i gh , Ml ow , U% I S2 , U2 ) 

dt/t 

1 f s(M '"s> h ' } M 2 t * duj/dt min(M 1 ,M 2 - Mi) 
~ 2 Js(M low) M(27r)i/2(5i-52) 3 / 2 M~ 2 

The M2 dependence of this rate can be seen by assuming the 
limits Mi om and Mhigh to scale with M2 , and approximating 
the dependence of S on M as a power-law S oc M a , where 
a = — (n + 3)/3 for the usual definition of the density power 
spectrum index n. The rhs of the above equation then scales 
as M 2 ~ a/ ' 2 . When n lies in the expected range between —2 
and —1, the mass dependence is very weak, roughly ~ M°' 2 . 
Taking into account that in the simulation one cannot really 
take infinitesimal time intervals, this dependence on final 
halo mass may be further weakened by the exponential term 
in the expression for /. 

Finally, we can also use these formulae to estimate the 
ratio of the growth rate through major mergers to that 
through "smooth" accretion (here defined as Mi < IM2). It 
is not easy to obtain analytic expressions for this ratio but 
it can easily be computed from the above formulae in the 
power law approximation for S(M). Here we give in Table 1 
the relative fraction for several typical values of a. Roughly 
speaking, major mergers are predicted to contribute 25% to 
30% of the total mass accretion, quite consistent with our 
simulation results. 
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Table 1. Table l:Relative growth rates due to major mergers and to other accretion modes. 
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Figure 6. Predicted relation between t * dui/dt (proportional to 
the dimensionless merger rate) and 1 + 2 for the concordance 
ACDM model (solid curve) and for an Einstein de Sitter universe 
(dashed curve) according to extended Prcss-Schechter theory. At 
redshifts above about 0.5 the two quantities are very nearly pro- 
portional to each other in the ACDM case also, as shown by the 
dotted straight line. 
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